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ABSTRACT 

The strongest spectroscopic dust extinction feature in the Milky Way, the 
broad absorption bump at 2175 A, is generally believed to be caused by aromatic 
carbonaceous materials - very likely a mixture of Polycyclic Aromatic Hydro- 
carbon (PAH) molecules, the most abundant and widespread organic molecules 
in the Milky Way galaxy. In this paper we report identifications of this absorp- 
tion feature in three galaxies at 1.4 < z < 1.5 which produce intervening Mgll 
absorption toward quasars discovered by the Sloan Digital Sky Survey (SDSS). 
The observed spectra can be fit using Galactic-type extinction laws, character- 
ized by parameters [R v , E(B - V)] ~ [0.7, 0.14], [1.9, 0.13], and [5.5, 0.23], 
respectively, where Ry = Ay/E(B — V) is the total-to-selective extinction ratio, 
E(B — V) = Ab — Ay is the color-excess. These discoveries imply that the dust 
in these distant quasar absorption systems is similar in composition to that of 
Milky Way, but with a range of different grain size distributions. The presence 
of complex aromatic hydrocarbon molecules in such distant galaxies is important 
for both astrophysical and astrobiological investigations. 

Subject headings: quasars: general, absorption lines-dust, extinction-individual 
(SDSS J144612.98+035154.4; SDSS J145907.19+002401.2; SDSS J012147.73+002718.7) 



1. Introduction 

The space between the stars of the Galaxy and external galaxies is filled with gaseous 
ions, atoms, molecules and tiny dust grains. These interstellar grains, spanning a wide 
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range of sizes from a few angstroms to a few micrometers, are important for the evolution of 
galaxies, the formation of stars and planetary systems. So far, one of the most best-studied 
properties of interstellar dust is its obscuration of starlight. The size and composition of 
interstellar dust are mainly inferred from, respectively, the extinction curve and its spectral 
features. The strongest spectroscopic interstellar extinction feature in the Galaxy is the 
broad 2175 A absorption bump. This feature was first detected by the Aerobee rocket 
observations (Stecher 1965). In this work we report the detection of the 2175 A dust exinction 
bump from individual intervening absorption systems in the spectra of three Sloan Digital 
Sky Survey (SDSS; York et al. 2000) quasars. 

The review of Draine (2003) provides an excellent summary of the properties of the 
2175 A dust extinction bump and theoretical constraints on its carrier. This feature is seen 
in extinction curves along lines of sight in the Milky Way (MW), the Large Magellanic Cloud 
(LMC), and some regions of the Small Magellanic Cloud (SMC). Extinction curves in the 
SMC bar region lack the 2175 A feature. The central wavelength of the feature varies by 
only ±0.46% (2a) around 2175 A, while its FWHM varies by ±12% (2a) around 469 A. 
Given its substantial band strength, the carrier responsible for the feature must contain at 
least one of the most abundant elements C, O, Mg, Si or Fe. Although the exact carrier 
is unknown, Draine (2003) states that "it now seems likely that some form of graphitic 
carbon is responsible" , most likely polycyclic aromatic hydrocarbons (PAHs; Joblin, Leger, 
& Martin 1992; Li & Draine 2001; Draine 2003 and references within). 

Several previous detections have been reported of the observed 2175 A feature in distant 
galaxies using quasar spectra. Pitman, Clayton, & Gordon (2000) presented a good summary 
of the state of the field as of three years ago, and ruled out several reported detections, which 
will not be discussed here. 

The highest redshift 2175 A bump observation to date is that of Vernet et al. (2001), 
who claimed a probable detection in a composite spectrum of radio galaxies at z ~ 2.5, 
although they acknowledged that the detection is difficult to confirm due to the low signal- 
to- noise ratio (SNR) and blending with Fell emission at A > 2300 A. Additionally, the dust 
responsible for the feature does not appear to be located in the host galaxies proper, but 
rather in the narrow emisson-line regions of these radio galaxies where the physical conditions 
may be quite different, perhaps leading to different carriers for the bump. 

The only previous detection of this feature from an individual intervening absorption 
system was that of Cohen et al. (1999), who detected the 2175 A feature in a damped Lya 
absorber at redshift z = 0.524 toward the BL Lac object AO 0235±164 at z = 0.94. They 
found a lower dust-to-gas ratio than in the Galaxy and mentioned that the 2175 A feature 
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"suggests there are differences from the average Galactic curve." 1 

Such differences appear to be rather common when the 2175 A feature is detected 
with data of sufficient quality to model the extinction curve involved. Malhotra (1997) 
detected the 2175 A feature in the composite absorption spectrum of 96 intervening Mgll 
absorption systems at redshifts 0.2 < z < 2.2. The strength of the feature was roughly 
consistent with a standard Galactic dust-to-gas ratio. Falco et al. (1999) reported detections 
in several z < 1 galaxies responsible for the gravitational lensing of background quasars. 
Many of the estimated extinction curves do not match the standard Ry = 3.1 Galactic 
curve (e.g., Ry = 7.2 ± 0.1 was found for a z\ = 0.68 spiral galaxy). Toft, Hjorth, & Burud 
(2000) detected the 2175 A feature from the lensing galaxy of B 1152+199 at z = 0.44 
and fitted its extinction curve by a Galactic-type extinction law with 1.3 < Ry < 2.1 and 
0.9 < E(B-V) < 1.1. Motta et al. (2002) detected a strong 2175 A bump in a lensing galaxy 
at z — 0.83. Their data are well fitted by a standard Ry = 2.1 ± 0.9 Galactic extinction 
curve, leading them to suggest that the lens galaxy of SBS 0909+532 contains dust like that 
of the Galaxy. Wucknitz et al. (2003) may have marginally detected the 2175 A feature in 
a lensing galaxy (later identified with a damped Lja absorber) at z — 0.93. They found 
that using extinction curves with a significant 2175 A bump reproduces the data better 
than curves without this feature. Recently Munoz et al. (2004) reported that the dust in 
the z\ = 0.68 lens galaxy ofB 0218+357 shows the 2175 A bump but produces a very flat 
ultraviolet extinction curve with Ry = 12 + 2. 

Given the range of extinction curves found in these extragalactic systems, large and 
well-selected samples extending to higher redshifts are needed to characterize the diversity 
and evolution of dust properties in the early universe. As a step towards this goal, we 
present the 2175 A feature identified in three individual intervening absorption systems in 
the spectra of three SDSS quasars. 

2. Data 

One of the goals of the SDSS is to obtain spectra for ~10 5 quasars to i = 19.1 {i = 
20.2 for z > 3 candidates), in addition to the ~ 10 6 galaxies that comprise the bulk of 
the spectroscopic targets (Strauss et al. 2002; Blanton et al. 2003). From astrometrically 
calibrated drift-scanned imaging data (Gunn et al. 1998; Pier et al. 2003) on the SDSS 
ugriz system (Fukugita et al. 1996), quasar candidates are selected primarily using color 



1 The average Galactic curve is parameterized by Ry = 3.1, where Ry = Ay/E(B — V) is the ratio of 
total to selective extinction and E(B — V) = A(B) — A(V) is the color excess. 
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criteria designed to target objects whose broad-band colors are different from those of normal 
stars and galaxies (Richards et al. 2002). Due to these inclusive criteria, the selection of 
candidates using i band magnitudes rather than blue magnitudes which are more affected 
by absorption and reddening, and its area and depth, the SDSS can reveal quasars with 
previously rarely seen properties. Relevant instrumental and observational details of the 
SDSS Data Release One (hereafter DR1) can be found in Abazajian et al. (2003). The 
DR1 spectra distributed by the SDSS cover approximately 3800-9200 A and have been sky 
subtracted, wavelength calibrated to the heliocentric frame, corrected for telluric absorption 
and also for Galactic extinction. One of us (P. B. H.) carried out a visual inspection of 
the SDSS spectra of the ~ 6350 DR1 quasars with z > 1.6 (Schneider et al. 2003). The 
goal was to identify examples of unusual quasar subtypes (primarily broad absorption line 
quasars, nitrogen-strong quasars, and dust-reddened quasars) for later comparison with the 
results of automated searches. This inspection uncovered several quasars with a possible 
2175 A signature, the most convincing three of which are SDSS J145907.19+002401.2, SDSS 
J144612.98+035154.4 and SDSS J012147.73+002718.7 (hereafter simply SDSS J1459+0024, 
SDSS J1446+0351 and SDSS J0121+0027). 



3. Fitting and Results 

The spectral index a redward of the Lyman a emission line — defined by fitting a 
power-law fx oc \~( a+2 } to the continuum — confirms possible reddening in these quasar 
spectra when compared to the composite quasar spectrum generated using a homogeneous 
data set of more than 2200 SDSS quasar spectra (Vanden Berk et al. 2001). As shown in 
Figure 1, compared to that composite's mean spectral index of « —0.46, SDSS J1446+0351 
has a « -1.93, SDSS J1459+0024 has a « -1.91, and SDSS J0121+0027 has a « -1.27, 
all suggesting either unusally red intrinsic continua or heavily dust-reddened spectra. We 
synthesized g, r, and % magnitudes from each spectrum independently and compared these 
values to the SDSS photometric measurements; this comparison demonstrated that the spec- 
trophotometric calibrations of the spectra are accurate. 

Since quasars have a range of intrinsic continuum slopes, to investigate the dust redden- 
ing properties in these absorption systems we simulate the dust-reddened quasar spectrum 
using template quasar spectra. The templates include the average SDSS quasar compos- 
ite spectrum from Vanden Berk et al. (2001), and the reddest-quartile and bluest-quartile 
composites from Richards et al. (2003) (composite spectra of SDSS quasars binned by con- 
tinuum color). We adopt different types of extinction laws and treat E(B — V), the amount 
of extinction, as a free parameter. We adopt the E(X — V)/E(B — V) relations from Pei 
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(1992) to calculate empirical extinction curves for the Milky Way, LMC and SMC. We also 
consider the "CCM Galactic extinction law" (Cardelli, Clayton & Mathis 1989), which gives 
a parameterized analytical form of A(X)/A(V) for wavelengths A > 1000 A: 

A(V) 1 ' + Rv 

where a(x) and b(x) are uniquely defined curves as a function of the wavenumber x = A -1 , 
and Ry = A(V)/E(B — V). (Note that the Milky Way empirical extinction curve used in 
Pei (1992) is equivalent to the CCM extinction law with R v = 3.08, the mean value from 
measurements along different lines of sight.) These extinction curves are then applied to 
redden the template spectrum. 

The simulated spectra are normalized to the data using the mean flux density between 
7000-7500 A or 7500-8000 A in the observed frame, since the red end of each spectrum is 
least affected by dust extinction. The reduced x 2 is calculated for each fit using the entire 
spectrum except for the emission line wavelength regions as defined in Vanden Berk et al. 
(2001) and for the wavelengths of narrow absorption lines. 

We now discuss the fitting results for the three individual quasars. Table 1 lists the 
redshifts of the quasars and of all absorption lines identified in the intervening systems 
associated with the putative 2175 A bumps. The narrow absorption line parameters were 
measured using the IRAF 2 standard software package splot. We assume that all dust along 
the line of sight is located at these intervening redshifts. For clarity, all the fitting attempts 
with different quasar composite spectra are summarized in Table 2; all the 2175 A measure- 
ments and the metallicity estimates for three absorption systems are summarized in Table 
3. 



3.1. SDSS J1459+0024 

The quasar redshift is z = 3.0124 ± 0.0005, and a strong intervening absorption system 
is seen at z = 1.3946±0.0001 (all absorption redshifts are measured from the Mgll doublet). 
The broad absorption feature centered between the Lya 1215 and Civ 1550 emission lines is 
consistent with being the 2175 A extinction feature from dust in the intervening absorption 
system at z — 1.3946. The restframe width of the broad absorption feature is ~ 400±60 A 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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FWHM, and its central wavelength is 2230±20 A at a redshift of z = 1.395. The measure- 
ments of FWHM and central wavelength are obtained with a Drude-like profile fit of the 
observed absorption feature 3 (the same measurements for §3.2 and §3.3). The la error is 
estimated by including different regions of continuum bracketing the feature into the fit (i.e., 
fit for several wavelength ranges between 2050-2550 A, in the dust restframe.). In the Milky 
Way this feature is observed to have an almost constant central wavelength Ao = 2175 ±9 A. 
This central wavelength difference between the Milky Way and this absorption system likely 
arises from the uncertainty caused by the Nv 1240 and Si IV 1398 emission lines bracketing 
the absorption feature. We rule out identifying this feature as a Civ broad absorption line 
trough because there is no evidence for broad absorption in other lines commonly seen in 
BAL quasars (e.g., Nv). 

By varying E(B — V) with the Milky Way, LMC and SMC empirical curves, it is 
found that the best fit to both the broad absorption feature and the large scale reddening is 
obtained using a Milky Way-type extinction curve w ith E(B -V) = 0.2. The 2175 A bump 
in the simulated spectra reddened with the LMC or SMC curves is not sufficiently deep to 
match the bump in the observed spectrum, indicating the existence of Milky Way-type dust 
in this absorber. To find a better fit to the flux level and the overall shape of the broad 
absorption, we add in Ry as another free parameter by adopting the CCM extinction law. 
With the CCM extinction law, the best fit is obtained with the reddest composite, which 
gives E(B -V) = 0.13 ± 0.01, Ry = 1.9lg| and a minimum reduced \ 2 ~ 1-8. The la 
error bar is estimated by variation of x 2 - With the average and the bluest composite, we can 
still obtain a good overall fitting with slightly larger \ 2 (Table 2, Fig. 2). The fitting with 
composite spectra of different colors indicates 0.7 < Ry < 1.9; the rather small value of Ry 
compared to the average value of 3.1 in the Milky Way (Savage & Mathis 1979) indicates 
that the size distribution of the dust grains is skewed towards substantially finer particles. 4 

Note that SDSS J1459+0024 is the only one of our three quasars detected by 2 MASS, 
and its J — K color of 1.22 ±0.1 7 places it in the reddest quartile of z = 3 quasars (Barkhouse 
& Hall 2001). The dust in the intervening galaxy would not have much of an effect at the 
rest-frame wavelengths probed by the J and K bands, so the red J — K color of SDSS 
J1459+0024 is good evidence that its UV-optical continuum is intrinsically red. This is 



3 The CCM functional form used a Drudc-likc profile to represent the 2175 A extinction hump: ~ x 2 j\{x 2 — 
x o) 2 + "f 2x2 ]i x = x o = Ag 1 , where Ao is the central wavelength and 7 is the FWHM. xo and 7 are 
determined by fitting the observed spectrum (Fitzpatrick & Massa 1986; Cardelli, Clayton, & Mathis 1989). 

4 Although the exact Ry value is dependent on the adopted composite spectrum type (reddest, average, 
or bluest), it is somewhat secure to conclude that this absorber system favors smaller Ry (than that of the 
Milky Way) and the SMC-type extinction curve without the 2175 A hump is not able to give a good fit. 
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consistent with our result that the best fit is obtained with the reddest quartile composite. 

3.2. SDSS J1446+0351 

The quasar redshift is z = 1.9452 ± 0.0023, and a strong intervening absorption system 
is seen at z — 1.5115 ± 0.0002. Only a few absorption lines are identified in this absorption 
system, as listed in Table 1. The quasar spectrum shows a broad absorption feature centered 
between the Civ 1550 and Cn] 2326 emission lines, consistent with the 2175 A extinction 
feature from dust in the intervening absorption system at z — 1.51. The width of the broad 
absorption feature is ~ 450±55 A FWHM, and its central wavelength is 2200±15 A at a 
redshift of z — 1.51. This small wavelength shift is probably due to the uncertainty produced 
by the C ill] 1909 line superposed on the absorption feature. 

With E(B — V) = 0.28 we obtain the best fit of the broad absorption feature and the 
large scale reddening with the empirical Milky Way-type extinction curve. Again, the 2175 A 
absorption in the simulated spectra based on LMC and SMC curves is not deep enough to 
match the observations. With the CCM extinction law, we obtain the best fit using the 
reddest composite with E(B - V) = 0.14 ± 0.01, R v = 0.7±° o i and a minimum reduced 
X 2 ~ 2.6. Using the average and the bluest composite, we obtained almost identical curves 
to the best fit with reddest composite and 0.5 < R v < 0.7 (Table 2, Fig. 3). This implies 
our fitting is not sensitive to the intrinsic quasar spectrum and the unexpected small value 
of Ry is not an artifact. The dust grain size distribution in this system also appears to be 
skewed towards smaller particles. 

3.3. SDSS J0121+0027 

The quasar redshift is z = 2.2241 ± 0.0001, and a strong intervening absorption system 
is seen at z — 1.3880 ± 0.0001. Standard lines identified in this absorption system are also 
listed in Table 1. Two other absorption systems are also seen in the spectrum: Civ, Siiv 
and possibly weak Mgll absorption at z — 1.955±0.001, and Civ, probably Nv and possibly 
Lja absorption at z — 2.200 ± 0.005. Again the quasar spectrum shows a broad absorption 
feature, centered between the Civ 1550 A and Cm] 1909 A emission lines, consistent with 
being the 2175 A feature from dust in the z = 1.388 intervening absorption system. The 
central wavelength in the rest frame for the broad absorption feature is 2250±25 A and 
the restframe width is ~ 400±60 A FWHM at a redshift of z = 1.388. The Civ emission 
probably causes the apparent shift of the dust absorption feature toward longer wavelengths 
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than seen in the Milky Way. 

Using the Vanden Berk (2001) SDSS composite spectrum and simply varying E(B — 
V) with the empirical curves, one cannot achieve a satisfactory overall fit to this quasar's 
spectrum. With the bluest-quartile composite spectrum from Richards et al. (2003) and 
E(B — V) — 0.14, we obtain the best fit of the broad absorption feature and the large scale 
reddening using the Milky Way empirical extinction. With the LMC and SMC empirical 
extinction curves, the 2175 A absorption is not sufficient to match the observed amplitude, 
also indicating Galactic-type dust content in this absorber (Fig. 4). 

We applied the CCM Galactic extinction law to the same bluest-quartile composite. 
We obtain the best fit with E(B — V) = 0.23^oai, Rv — 5-5 ± 0.1 and a minimum reduced 
X 2 ~ 1-3 (Fig. 4). Using the reddest composite of Richards et al. (2003) leads to difficulty 
matching the flux in the far-UV region. We estimate a range of 5.5 < Ry < 6.0 from the 
fitting results using the average and bluest composite (Table 2). The large value of Ry for 
this system implies that large dust grains are more prevalent than in the Milky Way diffuse 
ISM. 



3.4. Discussion 

The best fit parameters of both of the first two systems (SDSS J1459+0024 and SDSS 
1446+0351) include surprisingly small Ry values, implying a relative overabundance of finer 
dust grains compared to Milky Way dust. 5 In the Milky Way, so far the most extreme extinc- 
tion curve is found for the line of sight toward the HD 210121 high latitude translucent cloud. 
The extinction curve, with Ry « 2.1, is characterized by an extremely steep far-UV rise and 
by a relatively weak and broader 2175 A hump (Welty & Fowler 1992; Larson, Whittet, & 
Hough 1996). Detailed extinction, polarization and IR emission modeling implies that this 
region is rich in finer dust grains (Li & Greenberg 1998; Larson et al. 2000; Weingartner 
& Draine 2001; Clayton et al. 2003). Detailed modeling of the extinction curves of these 
intervening systems in terms of the silicate/graphite- PAHs model (Li & Draine 2001) are 
in progress (J.Wang et al. 2004, in preparation). Our goal is to infer their dust size distri- 
bution and composition. However, there are other possible explanations. For example, if 
the quasars are intrinsically redder than indicated by our fits, to reproduce their observed 
far-UV spectra would require extinction curves which are less steep in the far ultraviolet (i.e., 
which have higher Ry). We note that there are also uncertainties in the 5-6 /zm -1 range of 



5 The previously reported minimum Ry = 1.5 is associated with an elliptical galaxy at z — 0.96 (Falco 
et al. 1999). 
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the CCM extinction curve. 

Nonetheless, since the physical conditions governing dust properties evolve with redshift 
(e.g., star formation rate, metallicity, etc.), it should not be surprising that the dust prop- 
erties also evolve with redshift. If the diversity of grain sizes is real and not due to biases in 
the extinction modeling (in both our work and that of other groups mentioned in §1), it is 
extremely interesting for understanding dust formation at high redshift. Generally speaking, 
preferential removal of small dust grains (e.g. by coagulation in dense molecular clouds) will 
result in a gray extinction law with large Ry, while in star forming galaxies finer grains will 
dominate the dust size distribution as a result of dust destruction by grain-grain collisions 
or grain shattering by shock waves. 

Lack of the detailed knowledge of the metallicity of these absorber systems, we use 
the Milky Way gas-to-dust ratios iV(H i+R 2 )/E(B - V) = 5.8 x 10 21 cm" 2 mag -1 (Bohlin, 
Savage, & Drake 1978) and N(Rl)/E(B - V) = 4.93 ± 0.28 x 10 21 cm" 2 mag" 1 (Diplas & 
Savage 1994) to estimate the neutral gas content of these absorption systems. From our 
best fit E(B - V) values we estimate JV(Hi)« 6.4 x 10 20 cm^ 2 for SDSS J1459+0024 and 
N(R i)« 6.9 x 10 20 cm" 2 for SDSS J1446+0351. The strength of the Mg II absorptions suggest 
that those absorbers may be damped Lja systems, which have a lower limit to their column 
densities of A r (Hl)=2 x 10 20 cm -2 . From the SDSS spectra of the two quasars (spectral 
resolution ~2000), the equivalent width ratio of Mgll AA2796,2803 indicates that both lines 
are either saturated or nearly so. Therefore we can constrain the range of column densities 
for Mgll: 7.7 x 10 13 < JV(Mgii)< 3.1 x 10 17 cm~ 2 for SDSS J1459+0024 and 1.7 x 10 14 < 
iV(Mgii) < 7.0 x 10 17 cm -2 for SDSS J1446+0351. The lower and upper ends of these ranges 
represent the linear and damping regimes of Mgll, respectively. A lower limit estimate of the 
Mgll abundance compared to solar abundance gives [Mg/H]=-2.7 for SDSS J1459+0024 
and [Mg/H]=-2.2 for SDSS J1446+0351. No Znll lines were found in either system, so we 
cannot estimate the dust depletion. Some recent surveys have discovered strong molecular 
hydrogen lines in damped Lja absorbers with strong dust depletion, suggesting that the 
formation of H2 on the dust grains is the dominant formation process (Ge & Bechtold 1997, 
1999; Ge, Bechtold & Kulkarni 2001). HST observations of these quasars will be extremely 
valuable to verify this hypothesis. 

For SDSS J0121+0027, using the Milky Way gas-to-dust ratio we estimate an average 
hydrogen column density iV(Hl)^ 1.1 x 10 21 cm~ 2 from our best fit E(B — V) = 0.23. Again, 
this suggests that the strong Mgll absorption system is a damped Lja absorption system. 
Based on the equivalent widths of the weak absorption lines of Zn II 2026 and Fe II 2587 in the 
z = 1.388 absorber, we estimate the column densities for Znll and Fell are 3.4 x 10 13 cm~ 2 
and 1.6 x 10 14 cm~ 2 , respectively, assuming they are both on the linear part of the curve 
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of growth. Since both Znll and Fell are the dominant species for these metals in DLAs, 
the estimated Zn and Fe abundances are [Zn/H] ~ —0.2 and [Fe/H] ~ —2.6, respectively. 
The relative depletion between Fe and Zn, [Fe/Zn] ~ —2.4, is very large, similar to heavily 
depleted diffuse clouds in the Milky Way, such as the ( Oph cloud ([Fe/Zn] ~ —1.6, Savage 
& Sembach 1996). This dust depletion is the largest among all of the high-redshift DLAs 
searched for dust and molecular absorption to date (Ge, Bechtold & Kulkarni 2001; Ledoux, 
Petitjean, & Srianand 2003). The large dust depletion provides additional evidence that the 
strong dust extinction feature in this system is real. 

Although the precise nature of the carrier of the 2175 A extinction feature remains 
unknown, it is generally accepted that this feature arises in some types of aromatic carbona- 
ceous materials (Draine 2003; Li & Greenberg 2003). A proposal involving a cosmic mixture 
of many individual PAH molecules, radicals, and ions is receiving increasing attention (Li 
& Draine 2001). PAHs, the most abundant and widespread organic molecules in the Milky 
Way Galaxy as well as in other nearby galaxies (Allamandola & Hudgins 2003), reveal their 
existence in interstellar space by emitting a distinctive set of broad spectral lines at 3.3, 
6.2, 7.7, 8.6 and 11.3 /im (Li & Draine 2001). PAHs are thought to form predominantly in 
the atmospheres of carbon stars (Latter 1991; Allain, Leach & Sedlmayr 1996); from there, 
stellar winds will deposit PAHs into the surrounding ISM. 6 As the single largest repository 
of organic material in our galaxy, PAHs play an important role in prebiotic chemistry which 
may ultimately lead to the development of organic life. Under dense molecular cloud con- 
ditions, energetic processing (starlight irradiation and/or cosmic ray bombardment) of ices 
containing PAHs produces aromatic alcohols, ethers, ketones, and quinones; these are es- 
sential for important processes in living systems (Bernstein et al. 1999, 2002, 2003). Such 
complex organic molecules can be delivered to the surfaces of planetary objects via dust par- 
ticles (Clemett et al. 1999) and carbonaceous meteorites (Zenobi, Philippoz, Zare, & Buseck 
1989) — the formation mechanism of complex organic molecules found in such meteorites is 
not known for certain, but the deuterium enrichment seen in many such molecules suggests 
an interstellar origin (Cronin et al. 1993). 

The presence of PAHs in distant galaxies may thus be of astrobiological interest. The 
galactic environments in the galaxies discussed in this paper have probably evolved to the 



6 Wc should note that the origin and evolution of interstellar PAHs are not very clear. Other suggested 
sources for interstellar PAHs include (1) shattering of carbonaceous interstellar dust or of photoprocessed 
interstellar dust organic mantles (Greenberg et al. 2000) by grain-grain collisions in interstellar shocks (Jones, 
Tielens, & Hollenbach 1996); (2) in-situ formation through ion-molecule reactions (Herbst 1991). In the early 
universe, massive rate of star formation results in a much denser UV starlight intensity and a much greater 
supernova explosion rate. Both effects would enhance the destruction of PAHs. We would admit that even 
the origin and survival of PAHs in the Milky Way galaxy are not well understood. See Li (2004) for details. 
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stage where the above-mentioned organic molecules, and possibly even molecules of ribonu- 
cleic acid and of various proteins, can form on the surfaces of dust consisting of PAHs 
(Bernstein et al. 2002; Munoz Caro et al. 2002). If this is the case, then the building blocks 
of life may have been present billions of years before the formation of the Earth. 

4. Conclusion 

We report three direct spectroscopic detections of the 2175 A dust absorption feature 
in quasar absorption systems at redshifts 1.4 < z < 1.5. These are the first detections 
of this feature in individual Mgll absorption systems. From fitting the composite SDSS 
quasar spectra reddened by a CCM Galactic extinction law to the observed spectra, we 
derive best-fit reddening parameters E(B - V) « 0.14, R v « 0.7 for SDSS J1446+0351, 
E(B - V) « 0.13, R v « 1.9 for SDSS J1459+0024 and E(B - V) « 0.23, R v « 5.5 
for SDSS J0121+0027. The various R v values in these systems compared to the average 
Galactic value of Ry = 3.1 indicates a wide range of dominant grain sizes among intervening 
absorption systems. We found that although the intrinsic slope of the quasar continuum in 
each system is unknown, we can still meaningfully constrain the Ry value in each system. If 
the presently favored PAH model for the 2175 A feature carrier is correct, we have detected 
complex organic molecules in the young universe, about 9 billion years ago. The presence of 
PAHs at such large lookback times may be of astrobio logical interest. 
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Table 1. Strong lines identified in dusty intervening absorption systems 



Arcst (A) 


Ion 


Aobs(A) 


EW Icst ± 1<T error(A) 


/ 


z 


SDSS J1459+0024 zq = 3.0124 


2852.9642 


Mgl 


6830.6 


0.50±0.04 


1.81 


1.3942 


2803.531 


Mgll 


6713.6 


1.67±0.04 


0.295 


1.3947 


2796.352 


Mgll 


6695.7 


1.96±0.04 


0.592 


1.3945 


2600.1729 


Fell 


6226.4 


0.92±0.08 


0.203 


1.3946 


2586.6500 


Fell 


6194.6 


0.58±0.17 


0.0573 


1.3947 


2382.7652 


Fell 


5705.6 


1.09±0.12 


0.328 


1.3946 


2344.2139 


Fell 


5613.0 


0.50±0.04 


0.108 


1.3945 


1670.7874 


A1II 


4003.8 


4.97±0.88 


1.88 


1.3959 


SDSS J1446+0351 zq = 1.9452 


2803.531 


Mgll 


7022.5 


3.66±0.24 


0.295 


1.5113 


2796.352 


Mgll 


7041.8 


4.50±0.20 


0.592 


1.5117 


2600.1729 


Fell 


6532.5 


2.31±0.12 


0.203 


1.5120 


1670.7874 


A1II 


4196.03 


0.96±0.12 


1.88 


1.5114 


SDSS J0121+0027 z Q = 2.2241 


2852.9642 


Mgl 


6814.0 


0.92±0.08 


1.81 


1.3884 


2803.531 


Mgll 


6695.0 


2.72±0.08 


0.295 


1.3881 


2796.352 


Mgll 


6677.4 


3.52±0.08 


0.592 


1.3879 


2600.1729 


Fell 


6209.2 


1.84±0.13 


0.203 


1.3880 


2586.6500 


Fell 


6176.2 


0.54±0.04 


0.0573 


1.3877 


2344.2139 


Fell 


5597.7 


1.05±0.17 


0.108 


1.3879 


2026.137 


Znll 


4836.6 


0.50±0.04 


0.412 


1.3871 


1862.7895 


A1III 


4448.8 


0.50±0.08 


0.268 


1.3882 


1854.7164 


A1III 


4430.2 


1.05±0.13 


0.539 


1.3886 


1808.0126 


Sill 


4309.2 


0.75±0.08 


0.0055 


1.3834 


1670.7874 


A1II 


3989.6 


1.13±0.08 


1.88 


1.3879 



Note. — Vacuum wavelengths, oscillator strength / are from Morton (1991). Equivalent 
width is measured in the absorber restframe. 
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Table 2. Fitting log with CCM extinction curve using different composite spectra 



Composite used 


[R v , E(B-V)} 


Reduced \ 2 


SDSS J1459+0024 


Reddest 


[1.9, 0.13] 


1.8 


Average 


[0.7, 0.14] 


2.0 


Bluest 


[0.7, 0.14] 


2.0 


SDSS J1446+0351 


Reddest 


[0.7, 0.15] 


2.6 


Average 


[0.5, 0.16] 


3.0 


Bluest 


[0.6, 0.18] 


2.9 


SDSS J0121+0027 


K oddest 


[6.0, 0.17] 


5.8 


Average 


[6.0, 0.23] 


2.5 


Bluest 


[5.5, 0.23] 


1.3 



Note. — "Reddest" and "Bluest" refer to the reddest and 
bluest quartile composite quasar spectrum from Richards 
ct al. (2003). "Average" refers to the composite quasar 
spectrum using median combining from Vandcn Berk et al. 
(2001). 



Table 3. Summary of the 2175 A measurements and metallicity estimates for three 

sightlines 



Line-of-sight 


z abs 


A (A) a 


FWHM(A) b 


JV(Hl) (cm- 2 ) c 


[X/H] d 


SDSS J1459+0024 
SDSS J1446+0351 
SDSS J0121+0027 


1.3946 
1.5115 
1.3880 


2230±20 
2200±15 
2250±25 


400±60 
450±55 
400±60 


6.4 x 10 20 
6.9 x 10 2() 
1.1 x 10 21 


[Mg/H] ~ -2.7 
[Mg/H] ~ -2.2 
[Zn/H] ~ -0.2; [Fc/H] ~ -2.6 



a Central wavelength of the 2175 A feature. 

b Full Width Half Maximum of the 2175 A feature. 

c Estimated Hi column density in the absorption system along line-of-sight. 

d Estimatcd metallicity [X/H] = log(X/H)— log(X/H)Q in the absorption system. Reference solar abun- 
dances are from Savage & Sembach (1996). 
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1000 1500 2000 2500 3000 
Quasar Restframe Wavelength (A) 

Fig. 1. — Comparison between the composite quasar spectrum from Vanden Berk et al. 
(2001) (top; spectral index a = -0.46), the observed spectrum of SDSS J1459+0024 (second 
from top; spectral index a = —1.91), the observed spectrum of SDSS J1446+0351 (third from 
top; spectral index a = —1.93) and the observed spectrum of SDSS J0121+0027 (bottom; 
spectral index a = —1.27). Power law fits to estimated continuum flux are shown as the 
dashed lines. The broad absorption features at ~1400 A in the quasar SDSS J1459+0024 
restframe, at ~1900 A in the quasar SDSS J1446+0351 restframe and at ~1700 A in the 
quasar SDSS J0121+0027 restframe are identified as the 2175 A dust extinction feature. 
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Quasar Restframe Wavelength (A) 
895 1194 1492 1790 2089 2387 




h ... .I i .... i .... i .... i .... -I 
1500 2000 2500 5000 5500 4000 
Dust Restframe Wavelength (A) 



Fig. 2. — Comparison of the observed SDSS J1459+0024 spectrum (black) with model 
spectra obtained by reddening the SDSS composite quasar spectrum (top panel: the reddest 
quartile composite from Richards et al. (2003); middle panel: average composite from Vanden 
Berk et al. (2001); bottom panel: the bluest quartile composite from Richards et al. (2003)) 
with three types of extinction curves: Milky Way with various Ry (red), LMC (green), SMC 
(blue). All spectra are normalized so that they have equal flux density at 3200 A. The model 
spectra using the composite from Richards et al. (2003) can not cover the entire observed 
spectrum simply because the bluest end of the composite starts from ~ 1195 A. The best fit 
is obtained using the reddest quartile composite, reddened by a Milky Way-type extinction 
curve with Ry « 1.9. 
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Quasar Restframe Wavelength (A) 

1279 1705 2132 2558 2985 3411 




1500 2000 2500 3000 3500 4000 
Dust Restframe Wavelength (A) 

Fig. 3. — Comparison of the observed SDSS J1446+0351 spectrum (black) with model 
spectra same as Fig. 2 but for SDSS J1446+0351. All spectra are normalized so that 
they have equal flux density at 3400 A. The best fit is obtained using the reddest quartile 
composite, reddened by a Milky Way-type extinction curve with Ry ~ 0.7. 
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Fig. 4. — Comparison of the observed SDSS J0121+0027 spectrum (black) with model 
spectra same as Fig. 2 but for SDSS J0121+0027. All spectra are normalized so that they 
have equal flux density at 3300 A. The best fit is obtained using the bluest quartile composite, 
reddened by a Milky Way-type extinction curve with Ry « 5.5. 



